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1 .  Introduction 


Mixed  Poisson  processes  play  an  important  role  in  many  branches  of  applied 
probability,  for  instance  in  insurance  mathematics  and  physics  (see  Albrecht 
(1985)  and  Pfeifer  (1986)  for  recent  surveys).  They  belong  to  the  class  of 
elementary  pure  birth  processes  {N(t);t  >_ 0  )  with  standard  transition 
probabilities 


p  (s,t)  =  P(N(t)  =  m  |  N(s)  =  n),  0  £  n  <  m,  0  <_  s  <_  t,  (1) 


possessing  right- continuous  paths  and  positive  and  finite  birth  rates 


X  (t)  =  lim  i  p  .(t,t+h),  n,t> 

n  '  ,  h  *n,n+lv  — 

fitU 


o. 


(2) 


and  all  finite-dimensional  marginals  of  the  jump-time  sequence  { T^;  n  ^  0  } 
are  absolutely  continuous  with  respect  to  Lebesgue  measure  (see  Pfeifer  (1982)). 
For  such  processes,  the  jump  times  form  a  Markov  chain  with  transition 
probabilities 


P(T  >  t 
n 


n-1 


1  " 

=  s)  =  - 2 - , 

1-Fn(s) 


0  <  s  <  t,  n  >1 


and  initial  distribution  function  where 


1  -  Fn(x)  =  exp 


x  >  0,  n  >  0 


(3) 


(4) 


(the  F  are  in  fact  all  cumulative  distribution  functions),  hence  all  F 


are  absolutely  continuous  with  densities  f  (say),  and  the  conditional 
densities  for  the  transition  probabilities  can  be  represented  as 


f  (s) 

f  (s  |t)  =  - 2 -  ,  0<  t  <s,  n>  1. 

1 -*•„(*>  -  -  - 


Moreover,  the  birth  rates  coincide  with  the  hazard  rates 


f(t) 

Xn(t)  =  - - - 

1  -  Fn(t) 


a.e. ,  n,t  >  0. 


If  especially  {N(t);  t>  0}  is  a  mixed  Poisson  process,  then  also 


OO  .  oo 

An(t)  =  xn+1e  dG(x)yf  |  xni 


e  dG(x)  ,  n,t  >0, 


where  G  is  the  cdf  of  the  mixing  random  variable  A  (say).  In  fact,  Lundberg 
(1940)  has  proved  that  such  a  representation  characterizes  the  intensities 
of  a  mixed  Poisson  process. 

In  terms  of  random  variables,  a  mixed  Poisson  process  behaves  like  a  homogeneous 
Poisson  process  with  rate  A  given  A  =  A  ,  from  which  it  also  follows  that 


An(t)  =  E(  A  |  N(t)  =  n),  n,t  >  0. 


The  following  result  completes  some  of  Lundberg 's  (1940)  results  on  the 
asymptotic  behaviour  of  the  intensities  for  mixed  Poisson  processes. 


■-'.■vVvVv' 

••  .Vy-N. 


Lemma.  Let  { t^;  n  >  1  }  be  a  sequence  of  positive  real  numbers  converging 


to  t  >  0  such  that  |  ~  -  1  |  =  o(n” n  00  .  Then,  if  1/t  is  a  point 


of  increase  of  G,  we  have 


lim  X  (nt  )  =  — 
n  n'  t 
n 


3  2 

Proof.  Let  e  >  0  be  chosen  in  such  a  way  that  n  e  -►  0,  n  e  -* 
-  n  J  n  ’  n 

t 

(t—  -  l)/e  ->  0  for  n  00  .  From  relation  (7)  it  follows  that 

T>  II 

oo  g  °° 

X .  (nt  )  =  ~  en(xt  )n+1e_nxtn  dG(x)  f  en(xt  )ne-nxtn  dG(x) 
n  n  n  '  /  n 

11  A  /A 


(1+rn}/tn  /  (1+ 

~  J  xtnexp(-  J(l-xtn)2)  dG(x)/ 

11  ^n^n  '  (1~e- 


(1+e  )/t 

r  n"  n 


exp(  -  |(l-xtn)2)  dG{x ) 


(1-e  )/t 
v  n  n 


—  for  n  -+  °°  .  This  proves  the  Lemma, 

n 


The  above  result  is  in  general  not  true  without  further  conditions  on  G  as  can 
e.g.  be  seen  by  mixing  distributions  concentrated  in  a  single  point  X  >  0; 

here  X^t)  =  X  for  all  n  and  t,  X  being  the  only  point  of  increase  of  G. 

✓ 

As  an  example,  consider  a  Polya-Lundberg  process  where  A  follows  a  gamma 

.  2 
distribution  with  mean  y  >0  and  variance  aU  ,  a  >  0.  Here 


,  t  * \  _  1  +  an 

X  (t;  -  y  -  , 

1  +  ya  t 


n,t  >  0, 


from  which  the  validity  of  (9)  can  be  seen  explicitly,  for  all  t  >  0. 


2.  The  martingale  characterization 


Let  {'n(Tn_^);  l)  denote  the  sequence  of  post-jump  intensities. 

In  the  light  of  (2),  the  post-jump  intensities  describe  the  transition  behaviour 
of  the  process  immediately  after  a  jump  has  occured.  The  following  result 
gives  a  characterization  of  mixed  Poisson  processes  by  a  martingale  property 
of  this  sequence. 

Theorem  1  Let  {  N(t);  t  >_  0  }  be  an  elementary  pure  birth  process  with  intensities 
( \.(  t );  n,t  >0}  and  jump  times  (T  ;  n>  0)  .  Let  for  n>  1  denote  A  the 
cr-field  generated  by  Then  (N(t);  t>  0  }  is  a  mixed  Poisson  process 

iff  the  post-jump  intensities  {X  ^(T^  ^);  n  >1}  form  a  martingale  with  respect 
to  {A  ;  n  >  1 }  . 

Proof .  Due  to  the  Markov  structure  of  jump  times  the  martingale  property  of 
the  post-jump  intensities  is  equivalent  to 

E^An+l^Tn^  I  Tn-1  =  t)  =  a.s.  for  all  n  >1  ( 11 ) 


which  by  (5)  and  (6)  is  in  turn  equivalent  to 


Ws)  fn(s)  „  fn<‘> 

ds  =  -  a.e. 


t^W8)  X-Fn(t) 


saying  that  the  density  f  is  differentiable  a.e.  with 


f(t> .  - 

n  n 


f  ( t )  a.e, 


or  equivalently 


(12) 


(13) 


4  - 


w  W-r„,t  - -  -  =  af  log  f  (t)  a.e.  (U 

dt  n*l  l-FnU(t)  fn(t>  dt 

Integration  of  this  last  relation  shows  that  there  are  constants  c^>  0  such 


1-F  ^,(t)  =  c  f  (t),  t  >  0, 
n+1  n  n  — 


which  in  turn  implies  that  f  is  absolutely  continuous  and  the  recursive 
formula 


f  x,(t)  =  -  c  f'(t) 
n+1  n  n 


holds  everywhere  on  [O,00)  .  By  induction,  we  see  that  all  derivatives  of 
fR  exist  on  [  0,  °°)  ,  and  that  for  all  n  1, 


f  (t)  =  (-i)n  n  c  An'(t),  t  >  0. 
n  k=0  *  U 


Since  by  assumption,  the  intensities  (and  hence  all  f  )  are  positive  and 
finite,  we  have 


(-l)n  f^n)(t)  >  0,  n,t>0. 


The  density  f^  thus  is  completely  monotonic  on  [0,°°)  ,  hence  by  Bernstein's 
(1928)  theorem  there  is  a  bounded  and  non-decreasing  right-continuous  function 
H  such  that 


In  fact,  since  f  is  a  density,  we  have  that  ^  dH(x)  =  dG(x)  is  a  probability 
measure  from  which  it  follows  that 


f  (t)  = 
n 


n-1 

n  C1 
k=0  1 


f  n+1  -xt  \ 
x  e  dG( x ) 


1-F  (t)  = 
n 


n-1 
H  c 
k=0  1 


xne  dG(x),  n,t  >  0. 


(20) 


Hence  relation  (7)  is  satisfied,  saying  that  {  N(t);  t  >  0}  must  be  a  mixed 
Poisson  process  with  mixing  distribution  dG(x). 

Conversely,  since  every  mixed  Poisson  process  has  intensities  of  the  form 
(7),  it  is  easily  seen  that  relation  (12)  holds,  hence  the  post-jump 
intensities  possess  the  martingale  property,  which  proves  the  theorem. 


It  should  be  pointed  out  that  since  f^O)  <  00  by  our  assumptions,  the  mixing 
random  variable  must  be  integrable  with 


E(  A  )  =  f Q(0). 


(21) 


A  simple  application  of  the  Martingale  Convergence  Theorem  (see  e.g.  Billingsley 
(1979))  then  shows  that  the  post- jump  intensities  converge  a.s.  to  some 
integrable  random  variable  since  also 


fx(t)  fQ(t) 
1-FX( t ) 


dt 


fA(t) 

— - fn(t)  dt  =  f JO). 

f0(t)  0  0 


(22) 


The  question  now  is  what  the  possible  limits  of  the  post- jump  intensities  are. 
The  following  result  gives  an  answer  to  this. 


Theorem  2.  If  A  is  the  mixing  random  variable  of  the  process,  then  the 
post- jump  intensities  converge  a.s.  to  A  . 


Proof.  For  any  mixed  Poisson  process,  we  have  (n+l)/Tn  ■*  A  a.s.  by  the 
strong  law  of  large  numbers,  applied  to  the  Poisson  process  with  rate  A  , 
conditionally  on  A=A  ,  and  by  the  law  of  the  iterated  logarithm, 

|  X  -  1 1'  =  o(/  loglog  n)  =  o(n  )  a.s.  for  n  •*  00  .  Since  also  A 

is  a.s.  concentrated  on  the  points  of  increase  of  G,  the  cdf  of  A ,we  have 
by  the  above  Lemma 

A  n+i(Tn)  =  A  n+i((n+1)(Tn/n+1)))  ~(n+l)/Tn+A  a.s.,  (23) 

which  proves  the  theorem. 

For  instance,  if  A  is  concen orated  on  two  poiiio«  \^  <  with  mass  a  and 

1-a  each  (a  >  0),  then 


V2' 

t 

-1<vi*V 

A  n(nt)  -  \ 

av  x+(l-a)v2, 

t 

-  1(VV 

(24) 

■  V 

t 

l1(v1,v2) 

where  l(v1/v2)  =  (loSv2  "  log  V1^V2  ”V1^  aS  can  Seen  from  Lundber§ 
(1940),  relation  (10 8).  Since  always  l/\i>  <  1(v^a>2)<  1/ \^,  it  can  explicitly 
be  seen  that 


X  n+l^’V 


with  probability  a 
with  probability  1-a  , 


(25) 


-  7  - 


i-e-  *  n+l(Tn)  -  A  a‘s* 
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